Dimensional dependence of naked singularity 
formation in spherical gravitational collapse 



Roberto Giambo 

E-mail: roberto.giambo@unicain.it 

Dipartimento di Matematica e Informatica 
Universita di Camerino 
62032 Camerino (ITALY) 

Sara Quintavalle 

E-mail: sara . quintavalleOunicam . it 

International School of Advanced Studies 
Universita di Camerino 
62032 Camerino (ITALY) 

Abstract. The complete spectrum of the endstates - naked singularities, or 
blackholes - of gravitational collapse is analyzed for a wide class of iV-dimensional 
spacetimes in spherical symmetry, which includes and generalizes the dust solutions 
and the case of vanishing radial stresses. The final fate of the collapse is shown to be 
fully determined by the local behavior of a single scalar function and by the dimension 
N of the spacetime. In particular, the "critical" behavior of the = 4 spacetimes, 
where a sort of phase transition from black hole to naked singularity can occur, is still 
present if A = 5 but docs not occur if A > 5, independently from the initial data of 
the collapse. Physically, the results turn out to be related to the kinematical properties 
of the considered solutions. 
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1. Introduction 

Understanding singularities has always been one of the most intriguing issues in General 
Relativity since its beginning. The mathematical prediction that gravitational collapse 
may lead to singularity formation hugely increased the attention over the study of last 
stage of heavy stars' life. Problems related to strong density regions of spacetimes 
are in need of an ultimate answer already and, over all, a satisfactory formulation of 
Penrose's Cosmic Censorship Conjecture [^]: the causal character, and the endstate, 
of singularities arising from a dynamical process such as an indefinite collapse, is still 
one of the favorite test-bed for relativity. A great amount of work in this direction has 
been done in the case of spherically symmetric 4-dimensional spacetimes: a number 
of collapsing models have been analytically studied where, under suitable assumptions, 
the arising singularity is not completely hidden behind a horizon, also when the latter 
forms. For instance the pioneering work of Christodoulou showed that it suffices 
removing homogeneity assumption from the paradigm of gravitational collapse leading 
to black hole - i.e. Oppenheimer-Snyder solution. These cases of naked singularities 
have been intensively explored, in particular Tolman-Bondi-Lemaitre dust clouds (see 
and references therein), and vanishing radial stress models |T^. 



Recently, a class of new solutions have been found out 0, including the above 
as particular cases, where naked singularities generically appears as an outcome of 
collapse. Physically, they describe the gravitational collapse of a class of anisotropic 
elastic materials, and are characterized by a particular choice of the equation of state 
that, in a certain coordinate system, allows to reduce Einstein Field Equations to a 
quadrature. In this paper, we find a natural extension of this class of solutions to the 
case of general A^-dimensional gravitation theory. The importance of higher dimensional 
models goes up e.g. to Kaluza-Klein theories, superstring theory, and brane-world 
models - see in particular ^ , where a description of the world with more than four 
non compact dimensions is proposed. 

In this perspective, the present study is motivated by a number of earlier and more 
recent works on spherically symmetric higher dimensional spacetimes: 1^ extend 



earlier well-known results and properties of the four dimensional scalar field collapse; 



Vaidya-adS four dimensional solution is generalized in |18, B3] to higher dimensions 



adding extra gravity terms to the action functional. Far from being exhaustive, more 



references on the subject of higher dimensional collapse are p, ||, ||, |TT], |T^. In particular, 
the class of solutions that we find extend again vanishing radial stress models as dust 
1^, "M]. We will find the complete spectrum of endstates, analyzing if and how it is 



modified by the dimension of the spacetime N. In particular, naked singularities will be 
proved to survive in any larger dimension, despite earlier results contained in |T^, ^| - 
see discussion at the end. 

It is worth noticing that some criticism arose to singularities occurring in 
astrophysical sources modeled with continuous media in the past, due to the fact that 
one can construct situations in which Newtonian systems made out of continua develop 



Dimensional dependence of naked singularity formation 



3 



singularities. As a consequence, singularities in these models cannot be considered as an 
exclusive product of General Relativity. It is difficult, however, to assess to which extent 
this phenomenon denies validity to continuous models, although a simple remark once 
made by H. Seifert |]2^ may be of help: on taking this point of view, one could discard 



the big-bang of the standard model as being an artifact of Newtonian gravity, since 
Friedmann equation holds - formally unchanged - also for the Newtonian cosmological 
models. 

The paper is organized as follows: section || is devoted to derive and the class of 
exact solutions, and to illustrate briefiy some particular cases. Physical reasonability 
conditions will also be imposed to the solution, together with conditions that will ensure 
formation of singularities, whose endstate will be analyzed in section ^. In section ^ 
we will show how to complete the model, matching the solution to a suitable exterior 
spacetime. In the final section we discuss the results found, relating them to kinematical 
properties of the spacetime. 

2. The solution in area— radius coordinates 

The general spherically symmetric line element in comoving coordinates {t,r,6^), i = 
1, . . . , N — 2, is given by 

= -e^-ii'r) ^^2 ^ ^^^^ ^^-1 ^^2 ^ ^^^^ ^)2 ^^2^_^^ 

where dfi^„2 = E.=I^(n}=\ sin^ ^^■) (d^^^^ jhe source of the gravitational field will be 
given by an elastic material in isothermal conditions. Generalizing the = 4 case , the 
property of the source are encoded in a state function depending on the space-space part 
of the metric, that is - using spherical symmetry assumption - w = w{r, R, rf) [|^, [l^ . 
The stress energy tensor is given by 

^ = ® ^ + pAr ^-^+Pt {de^ ® ^) , (2) 

where, introducing the matter density 

p= {N ~2){8nE)-^^R-^ 

[E is an arbitrary function of r), the internal energy e and the stresses Pr and pt are 
given in terms of the state function by 

dw 1 dw , . 

Although the comoving coordinates usually yields the natural system to describe 
the physical evolution of the collapse, for our purposes, however, it will be convenient 
to introduce the area-radius coordinate system {r,R,6'^), first introduced by Ori EH] in 
the study of 4-dimensional charged dust, in such a way that (|I]) becomes 

ds^ = -Adr^ - 2B dRdr - u-^ dR^ + R^ dQ%_^, (4) 
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with A,B,u unknown functions of {r,R). In this way the internal energy will depend 
only on one field variable, t], and on the two coordinates r, R. We introduce the function 

A = - Au-^ = 

so that Einstein field equations G^. = 8nT^, = SnT^ and = 8nT^ can be 
expressed in terms of A, A and u as follows: 

(1 - ^)[{N - 3)(1 - A/ A) - R{AIA)^n\ = 8vri?V, (5) 



;i - —)R-\A/A),r = -SvrVA + Au-^u'\e + p.), (6) 



/u^ + A/A{VA),n-{u-'),. = 0. (7) 

Equation can be integrated to give A/ A in terms of Pr{r,R,ri). Therefore, if one 
removes dependency on the comoving field variables, assuming that Pr in satisfies 

dprjr, R,ri) ^ ^ 
df) ' 

or equivalently 

w = h{r, R) + i{r, R)7]-^/^ (8) 
with h, i arbitrary, then one obtains 

A = A{1 - 2^ R^-^) (9) 
where, in view of (|]) and (P), \1/ is the function 

^(r, R) = F{r) + T i{r, a) da, (10) 

with F{r) arbitrary function of r, and Ro{r) describing r at initial (comoving) time, 
that will be chosen equal to r hereafter. The function \1/ ([T0|) is Misner-Sharp mass of 
the system, defined by the relation 1 — 2\1'_R^~^ = g(VR, V-R). Now, inroducing 

y(r, R) = E{r) ^„ (r, R) h{r, i?)"\ (11) 
the field equations (||)-(0), in view of (||), (|) and (p!0D , simply become respectively 

^2 = 2^/?3-A^ - 1 + (12) 

and {\^)^R + Y-\u-^),r = 0, that can be integrated, using the initial condition, to give 

^(r, R) = r da + ] (13) 

Jr Y {r,a)U[r,a)'^ Y[r,r)u[r,r) 

Then, we conclude that the class of exact solutions found expresses all the metric 
unknown functions in (^ in terms of two arbitrary functions (\1/, Y) of r and R. 

We stress the fact that the constitutive function w{r,R,ri) as equation of state, 
introduced at the beginning of this section, uniquely and completely carries on the 
physical properties of the matter, regardless of possible anisotropics. Isotropy of the 
matter is characterized when w can be written as a function of the matter density p 
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only. In this case, = pt and both can be seen as a function of the energy density e only, 
as one can easily calculate from (|^). When w fails to be a function of p only, anisotropy 
comes into play, but it is not needed any other relation to close the system, because 
of equations (^. Another way to see this is to observe that equations (3) identically 
imply the conservation law arising from one of the Bianchi identities written in comoving 
coordinates, and again this is of course an outcome of having assumed that the source 
is an elastic continuum in isothermal conditions. Of course, the requirement given by 
(1) is exactly the state function characterizing the class of function considered, and the 
fact that the arbitrary functions can be viewed in terms of the kinematical properties of 
the continuum is a very well known consequence of the structure of the field equations 
within the assumed symmetries and holds for all the models of this kind. 



2.1. Examples 



The components of the stress energy tensor are generically nonzero, as readily calculated 
from (P|)-(pi]D, and are given by 



Y 



Pr 
Pt 



87ri?^-2 ' 



OR ^ ' 



From these expressions we can recognize some particular cases: 
(i) dust spacetimes |T^, occurring when both and Y are functions of r only; 



(14) 
(15) 
(16) 



(ii) vanishing radial stress solutions pO|, that occur when \E' is a function of r only, but 
Y may also depend on R, 

(iii) acceleration free solutions, when y is a function of r only but \E' may also depend 
on R (note that the norm of the acceleration is simply given by F^j) 



2.2. Energy condition and shell focussing singularity occurrence 

On the above class of solutions, some conditions will be imposed, as requirements on 
and Y. Since we want to obtain global gravitational collapsing models, we will consider 
the interior metric @ as defined on a right neighborhood [0, rt] of r = 0, for some 
Tfo > 0, and match the above solutions at r = r?, with some exterior spacetime to be 
defined later (see section For this reason, in the following we will consider \E' and Y 
as defined on the set {{r,R) : r G [0,rb], R G [0,r]}. 

As a physical reasonability condition, WEC on the metric (^ will be required, but 
in view of (p^)-(p!6|), it suffices that 

{N - 2)^„ Y-' > R{^,r Y-'),R , {N - 2)^,n > R^j^n 
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Moreover, we impose the condition of decreasing initial energy, i.e. eo{r) := e(r, r) must 
be a decreasing function of r: 

(A^ - 2) 

^,^^(r, r) + ^,R/j(r, r) + 2^,,K(r, r) < [^,r(r, r) + ^^^(r, r)]. 

The functions \1/ and y must be chosen in such a way that the spacetime is regular at 
initial (comoving) time, and a (shell focussing) singularity forms, for each shell r G [0, r^], 
in a finite amount of time. Therefore, first of all shell crossing singularity formation must 
be avoided, and to this aim it must be required that > when R> 0. By inspection 
of (|T3D, sufficient conditions for this to happen are given by 

VA(r,0)>0, u^r{r,R)>0, Vr g]0, r^], i? g]0, r]. 

Moreover, it must be observed that the use of the (r, R, 6'*) coordinate system has 
the obvious advantage to parameterize the singularity with the straight line R = 0, but 
the drawback that both the regular and the singular centre are mapped into the point 
r = R = 0, and then it does not make a distinction between them, unless one does 
not consider the inverse function t = t{r, R). The function R, the derivative of R w.r.t 
comoving time, satisfies the identity u = — i?e~^, that can be formally integrated to give 
t{r, R) = J^e~'^u{r, a) da. Although the integrand yet contains an unknown function in 
the comoving coordinates, a key remark at this stage is to observe that e'" is bounded 
in a neighborhood of the centre, which allows to express the above conditions in terms 
of \1/ and Y: it suffices that the function R^'^v^ is Taylor expandable at the centre 
(r = = 0), with expression given by 

i?^-V= hi/R^ + E hi/W + ... (17) 

iJ^j=N~l i+j=N-l+p 

In particular, for the centre to become singular in a finite time, it must be required that 
{hN-ifl, • • • , 7^ 0. Hereafter, we will suppose, as already done in 

a := hN-ifi 7^ 0. 

Although this is a generic assumption, the results we are going to state can also be 



extended to the degenerate case a = 0, as done in 29 1 for A^ = 4 



3. Naked singularity vs. black hole formation 

The endstate of the singularity for these models will be studied. First, let us observe 
that the central singularity is the only one that can be naked. Indeed, under the above 

1 iV-l N-1 

assumptions, the apparent horizon Rh{r) is such that Rh{r) = a'^-^r^-^ +o{rN-3), and 
moreover, if th{r) and ts{r) are the comoving times when the shell labeled r becomes 
trapped and singular, respectively, then limr^o+ th{r) — ts(r) =0. 

To analyze the endstate of the central singularity we will study existence of null 
radial geodesies Rg{r) emanating from the (singular) centre, such that Rg{r) > Rh{f) in 
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a right neighborhood of r = 0. To do that, we will use a remarkable property of Rh{r), 
to be a supersolution of null radial geodesic equation 

d /? 

^ = uVA(Y -u). (18) 
dr 

Therefore, to have existence of such a Rg{r) as above, we will actually look for 

N — l 1 

suhsolutions of (|T8|) of the form RAr) = xr 'v-s with x > a'v-s that therefore emanate 



from the singular centre - so that Rg also will. Incidentally, this also explains why it 
suffices to look for radial curves: indeed, the projection of a nonradial geodesic on the 
(r, R) plane would be a supersolution of ([T^), so if the singularity is nonradially naked, 
is also radially naked. 

As it happens for the = 4 case, the endstate of the singularity is related to Taylor 
expansion of the function 

v^(r,0) = er"-' + o(r"-^), (19) 

but also the dimension N of the spacetime will play now a crucial role. Indeed, the 
condition for the existence of R^ as above is equivalent to the existence of x > a^f^ 
satisfying 

N — l 2 , I a , , I a „ 1 N{N-3) ^ 

-xr—^ < (1 - ./-_)(. /-_er"-i + xr^). (20) 



The above inequality gives the complete spectrum of the endstates since it provides a 
necessary and sufficient condition for the singularity to be naked. Indeed, if = 4 one 
recovers the well known results of that the inequality holds - and hence the singularity 
is naked - if n = 1,2, and if n = 3 a critical case happens when the endstate is related 
to the value of ^ in ([T9|), since it must be 2^ > (26 + 5\/3)a for the singularity to be 
naked. In larger dimensions, the singularity is naked if n = 1, VA^, and if n = 2, N = 5, 
provided 2^ > 21 ^fa. In all other cases a black hole forms. Then we observe that the 
critical behavior, when a phase transition from black hole to naked singularity occurs, 
depending on the value of ^, is a feature of dimensions A^ = 4, and A^ = 5, and is 
forbidden at larger dimensions. As one can see, the contribution of the dimension A^, 
when it is larger than four, basically enters in the behavior of the apparent horizon, 
that behaves like 7-1+2/(^^-3)^ which is no more an integer power of r as A^ > 6, and it 
is always leading upon the "kinematical" contribution of A^ - i.e. the last term in (|20|) . 
Since the contribution of V^{r, 0) is always an integer power of r - see below - this fact 
results in the lack of critical case when A^ > 6. 



4. Exterior spacetime and matching conditions 

In this section we will see how to complete the model, matching the interior solution 
studied so far with an exterior spacetime, and requiring that Israel-Darmois junction 
conditions hold along the matching hypersurface S = {r = rj,}. From (^) we observe 
that radial pressure Pr does not vanish in general along S, so we cannot expect to match 
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Figure 1. Asymptotic behavior of the apparent horizon with respect to the dimension 
near the centre. In the (r, R) plane, the shaded region represents the evolution of the 
solution. Fixed the value for a, the higher N is, the bigger is the trapped region lying 
between the r-axes and the horizon. 



the solution with a Schwarzschild exterior. In this case a natural choice for the exterior 



metric can be given by generalized Vaidya solutions ||Tl|, that for generic read 

2M{V,Ry 



ds: 



ext 



- 1 



dV -2dRdV + dn 



2 

N-2i 



and Israel-Darmois junction conditions simply become requirements on the mass 
function M(V, R) on the junction hypersurface. To find the conditions, it is 
convenient to work with the general interior metric written in co moving coordinates (|lD . 
Parameterized E with coordinates (r, 6'*) ^ (r, r^, 6'*), The first and second fundamental 
forms of S w.r.t. this metric read 



lL = -e2'^dr2 + i?Mn^_2, 



II: 



int 
S 



R i?'df2^_2 I , 



(21) 
(22) 



where a dash and a dot denote derivatives w.r.t. r and t respectively, and all functions 
are intended evaluated in {T,rb). Injection of S into the exterior spacetime reads 
in coordinates as {V{t),Y{t),9''), where l^(r),F(r) must be determined. The first 
fundamental form of S takes the form 

■ 2M{V{T),Y{r)) 



I 



V{Ty + 2V{t)Y{t) 



dr 



2 ' Y{TfdQ%_,, (23) 



r(r)^-3 

where, with a slight abuse of notation, we denote by a dot the derivative w.r.t. r. 
Comparing ( PTf ) with (|2^) gives 

Y{r) = R{T,rb), (24) 
2M 

Using these relations we can express the second fundamental form of S w.r.t. the exterior 
metric as 

1 



+ 2VY = e 



2u 



(25) 



11, 



ext 



-y(v- -x,yV' 



+ V[Y + x,yYV + XxV' + \xyV^ 



dr^ 
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-YiY + xV)dn%^,}, (26) 

where x = 1 ~ 2MY^~^ . Comparing angular terms in ( p2D and (^) and using (p5|) 
gives 

X = R"v - (i?e-'^)^ (27) 

that is continuity of Misner-Sharp mass across S. Therefore we find the differential 
equation for V{t): 

V{t) = —. . (28) 

R'r] + Re-^ 

At this stage, it remains to compare dr^ terms in the second fundamental forms. But, 
with some algebra, the above relations together with field equation R! = XR' + v'R 
simply reduce the condition to 

My(\/(r),r(r)) = 0. (29) 

Therefore, we conclude that generalized Vaidya solutions can always be matched to a 
spherically symmetric interior metric (|I|) along S, provided that conditions (0) and 
(p8|) hold, and the mass function MiV^Y) satisfies (pT]) and (p9|) on E. 

The above fact can easily be translated in area-radius formalism: it suffices to 
parameterize E with coordinates (cr, 6'*) ^ {rb,a,6^). In this case the injection of S in 
the exterior spacetime reads (V^(cr), cr, 6**), where V{a), in view of (p^-(p9D, becomes 

dV 1 

da u{rh,a){u{rh,cr) -Y{rh,a))' 
and the mass function satisfies 



(30) 



M(\/(a), a) = ^(r,, a), My(y(a), a) = ^,,j(r,, a). (31) 

It can be observed that an interesting subclass of the above exterior metric is given 
by the anisotropic generalizations of deSitter spacetime which is obtained taking 
M = M{Y). Obviously, in this case condition (^9]) is trivially satisfied, and (|3T|) simply 



reduces to require continuity of the mass across the junction hypersurface (see also 0). 
5. Discussion and conclusions 

There have been previous works trying to explain the endstate in terms of the 



kinematical properties of the spacetime, in particular the shear at initial time |Ty, [20[] . 
In the following we are going to address this point, relating the indices n,^ coming from 
(p!9|) to all kinematical properties (see also 0). The function \/A{r,0) can be split 



in the sum of Ii{r) + hi^), where 



/i(r) := ^ [ / . da, 

Y (r, r) or Jo u[r, a) 

^'^"■^ lo [ym ~ Y^) {'^Im) 
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The behavior of these quantities near the centre can be studied, to find that Ii{r) = 
parP~^ + o(r^~^), where p is given in ([T^), and a G R depends on the coefficients hij of 
order — 1 and N — 1 + p. 

Introduced the polynomials -Pfc('r) = Yl!j=Qhk^jjT\ the value of a is given by 

P^-i+p(r)rV(^-3) ^ 

On the other side, /2('") = ^^'^ + o(r^), where 6 G R and q is the order of the first 
nonvanishing term of (r, R) expansion at the centre. Then, n in ( jlQ]) is given by the 
smallest between p and q+ 1. Now, the shear of the solution can be controlled by the 
scalar 

<y = ^<y^''<y,u = -u^J{N/2 - 1)/{N - l){\og{R{u^^r')),n 
and on the initial slice R = r behaves like pa^r^ + ©(r^), where ctq = 



(Ar/2- l)/(iV- l)(P/v-i+p(l))/(4Piv-i(l))^/^- We deduce that the asymptotic 
behavior of the shear near the regular centre can be responsible for the quantity 
/i(r) only, and does not even control the value of the parameter a in the critical case 
(n = 2, = 5) - not to tell that one can conceive cases when ao = 0. On the other side, 
the norm of the acceleration is given simply by Y,fi, and then it rules the quantity l2{f)i 
but again the knowledge of the initial acceleration could not be enough to establish the 
value of q. We can conclude that the evolutions of both acceleration and shear infiuence 
the endstate of the gravitational collapse, but none of them can be considered as a 
stand-alone responsible, as the function \/~K{r, 0) is, together with the dimension A^. 

We observe that, if > 6, the singularity is naked only when n = 1. This 
is not in contrast with [10, 20 1, where dust and vanishing radial stress solutions are 



shown to produce a black hole when A^ > 6. Indeed, the special case considered in 
those paper are acceleration free, or more generally such that -\/A(r, 0) behaves like 
/i(r) anyway, and so the endstate is related to the first nonvanishing power of 



after the (A^ — l)-th order. In the cases produced in [0, |2^ the expansion for both 
\l/(r)/r^~^ and Y"^ is assumed to contain only even order terms, which excludes the 
possibility n = 1. Instead, in the case studied in the present paper a more general 
situation is considered, when may contain both odd and even order terms, but 

only terms of order A^ — 1 + 2k, A; G IN, when restricted on the initial slice R = r. 
In other words, solutions may be produced, when vl>/r^~^ and Y"^ are even at initial 
time, but later they evolve to allow also for odd order terms. All in all, the conclusion 
stated for A^ = 4 in is confirmed at higher dimensions, that the formation of naked 
singularities or black holes weakly depends on the initial data, but is essentially a local 
phenomenon, depending on the Taylor expansion of a kinematical invariant near the 
centre. The contribution of the dimension is basically related to the behavior of the 
apparent horizon, that forbids occurrence of critical cases when A^ > 6, and restricts, 
but still allows for naked singularity formation at any dimension. 
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